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Abstract 

Background: Susceptibility to Fusarium wilt disease varies among wild accessions of Arabidopsis thaliana. Six 
RESISTANCE TO TUSARIUM OXYSPORUM (RFO) quantitative trait loci (QTLs) controlling the resistance of accession 
Columbia-0 (Col-O) and susceptibility of Taynuilt-0 to Fusarium oxysporum forma specialis matthioli (FOM) are 
detected in a recombinant population derived from a single backcross of the F, hybrid (BQ)- In particular, the RF01 
QTL appears to interact with three other loci, RF02, RF04 and RF06, and is attributed to the gene At1g79670. 

Results: When resistance to FOM was mapped in a new BQ population, in which the loss-of-function mutant of 
At1g79670 replaced wild type as the Col-0 parent, RFOVs major effect and RFOTs interaction with RF02, RF04 and 
RF06 were absent, showing that At1g79670 alone accounts for the RF01 QTL. Resistance of two QTLs, RF03 and 
RF05, was independent of RFO 1 and was reproduced in the new BQ population. In analysis of a third BC q popula- 
tion, resistance to a second pathogen, F. oxysporum forma specialis conglutinans race 1 (FOCI), was mapped and 
the major effect locus RF07 was identified. 

Conclusions: Natural quantitative resistance to F. oxysporum is largely specific to the infecting forma specialis 
because different RFO loci were responsible for resistance to FOM and FOCI . The mapping of quantitative disease 
resistance traits in BQ populations, generated from crosses between sequenced Arabidopsis accessions, can be a 
routine procedure when genome-wide genotyping is efficient, economical and accessible. 

Keywords: Arabidopsis thaliana, Fusarium oxysporum, Quantitative disease resistance, QTL mapping, Crossover 
interference 



Background 

Fusarium wilt of Arabidopsis thaliana is an ideal patho- 
system for mapping, identifying and characterizing genes 
responsible for host resistance to vascular wilt fungi. A. 
thaliana, which is the preeminent subject of plant mo- 
lecular genetic and genomic studies, is susceptible to in- 
fection by three phylogenetically-distinct pathogenic 
forms, or formae speciales, of the soil-borne fungus Fu- 
sarium oxysporum [1,2]. In the field, F. oxysporum forma 
specialis conglutinans (FOC), F. oxysporum forma spe- 
cialis raphani and F. oxysporum forma specialis mat- 
thioli (FOM) are isolated from diseased Brassica species, 
radish (Raphanus sativus) and garden stock (Matthioli 
incana), respectively [3]. Fusarium wilt of A. thaliana 
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recapitulates the development of disease symptoms in 
field hosts [1]. 

The response of different accessions of A. thaliana to 
different formae speciales varies from complete resistance 
to ready susceptibility [1]. For example, the standard la- 
boratory accession Columbia-0 (Col-0) is completely re- 
sistant to FOM but expresses only partial resistance to 
FOCI. Taynuilt-0 (Ty-0), on the other hand, is susceptible 
to FOM but also expresses partial resistance to FOCI. 

Two strategies are used to map genes responsible for 
phenotypic variation in populations [4-6]. When the popu- 
lation of interest is wild and results from an indeterminate 
number of undefined crosses, a genome-wide association 
(GWA) study uses evidence of linkage disequilibrium to 
associate sequence polymorphisms within or near the 
genes responsible for the trait. Enabling GWA studies in 
the plant A. thaliana is the primary motivation for the 
1001 Genomes Project, which has generated whole gen- 
ome sequence for hundreds of wild accessions of A. 
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thaliana [7,8]. Indeed, the detection of functional sequence 
diversity in A. thaliana using GWA is reported [9,10]. 
However, GWA studies rarely detect more than a modest 
fraction of the sequence diversity responsible for variation 
in existing populations of plant and animal species [5,9,11]. 

Genetic linkage may be used to map the genes associ- 
ated with a trait to chromosomal intervals. However, this 
approach requires that the studied population is derived 
from controlled crosses between defined parents; and, 
only the genetic diversity distinguishing the parents of 
crosses is detected. Nevertheless, linkage analysis has 
been a powerful and successful approach for detecting 
and defining the genes responsible for complex traits in 
A. thaliana [12]. 

With plant species that readily inbreed, such as A. 
thaliana, recombinant inbred (RI) populations are al- 
most exclusively used to map and define genetic loci 
underlying natural traits [12,13]. RI populations in their 
simplest form originate from an outcross between parents 
with dissimilar genotypes. Unique recombinant genotypes 
of the parents are captured in dozens to hundreds of RI 
lines that result from propagating individual F 2 offspring 
by self-fertilization and single-seed descent. After several 
filial generations of inbreeding, RI progeny become largely 
homozygous and thus true-breeding RI lines. However, 
the effort to propagate and curate an RI population with- 
out introducing selection represents a substantial invest- 
ment in time and effort before QTL analysis begins. The 
effort to generate an RI population is offset by the fact that 
RI lines are immortal and can be retested innumerable 
times and reused in separate studies but need to be geno- 
typed just once. There are now dozens of published RI 
populations from crosses between wild accessions of A. 
thaliana [12,14-16]. Recently, a technique for generating 
haploid A. thaliana has made the generation of doubled 
haploid lines possible [17]. Like RI lines, doubled haploids 
are homozygous and thus immortal but require fewer gen- 
erations to create. 

Other mating strategies generate recombinant map- 
ping populations in less time and with less effort than it 
takes to generate RI lines. In particular, BCi populations 
are generated from crosses in two successive generations. 
An initial outcross between parental genotypes produces 
the Fx hybrid, which is then backcrossed to its recurrent 
parent. Each resulting BCi hybrid inherits a set of non- 
recombinant chromosomes from the recurrent parent and 
a set of recombinant chromosomes from the Fx hybrid. 
Because crossovers resulting from single meioses can be 
unambiguously assigned to recombinant chromosomes, 
the BCi mating scheme is often used to generate a model 
population for the evaluation of novel approaches to QTL 
analysis [18-20]. In addition, backcrossing is a common 
feature in traditional breeding schemes that seek to intro- 
gress new traits into elite crop varieties [21]. 



The appeal of BCi populations is undermined by the 
need for extensive genotyping, and very few studies of 
natural traits in A. thaliana have used BCi populations 
for genome-wide mapping [1,12,22]. Because each BCi 
hybrid possesses a unique recombinant genotype, it is 
necessary to genotype each tested BCi hybrid genome- 
wide. Without whole genome sequence information for 
the parents of a BCi population, the discovery of sequence 
polymorphism and their development into an appropriate 
set of DNA markers for genome-wide mapping is a time- 
consuming and laborious process. 

Nevertheless, prior genetic analysis of a BCi popula- 
tion shows that the qualitative resistance of Col-0 to 
FOM is a polygenic trait [1]. Six RFO QTLs, accounting 
for the resistance of Col-0 and susceptibility of Ty-0, 
segregate in a population generated by crossing Col-0 
and Ty-0 and then backcrossing the resistant Fi hybrid 
to its susceptible parent Ty-0. Among RFO loci, RFOl 
has the strongest association with resistance to FOM. 
RFOl also appears to interact with three other RFO loci, 
namely RF02, RF04 and RF06, because the three inter- 
acting loci have significant association only when recom- 
binant BCi hybrids also inherit the Col-0 allele of RFOl 
(RFOl-Q. RF02 is a receptor-like protein (RLP) gene 
that is homologous to the PSY1 peptide receptor gene, 
PSY1R [23]. The RFOl-lmked gene Atlg79670 is named 
RFOl because the Col-0 sequence of Atlg79670, as a 
transgene, enhances the resistance of Ty-0, and the loss- 
of-function allele of Atlg79670 (rfol) compromises the 
resistance of Col-0 [1]. Atlg79670 is a member of the 
wall-associated kinase-like kinase subfamily of receptor- 
like kinase (RLK) genes. 

Here, I map Fusarium wilt resistance in two new BCi 
populations (i) to address whether Atlg79670 alone is 
responsible for resistance attributed to the RFOl QTL, 
including interactions with RF02, RF04 and RF06, and 
(ii) to examine whether the same or different RFO QTLs 
mediate resistance to different formae speciales of F. oxy- 
sporum. In doing so, I present a methodology for genome- 
wide genotyping that makes the mapping of complex 
quantitative traits a routine procedure. Importantly, be- 
cause whole genome sequence is now available for most 
studied accessions, the same approach could be applied to 
crosses between any pair of Arabidopsis accessions. 

Results 

Resistance to FOM in rfol 

In prior mapping of resistance to FOM, RFOl was the most 
significant of six RFO loci in A. thaliana, and RFOl was 
epistatic to, or enhanced the resistance of, three other RFO 
loci [1]. In theory, the RFOl QTL could represent one gene 
or multiple genes. To appreciate whether Atlg79670 is re- 
sponsible for all or part of the resistance attributed to the 
RFOl QTL, resistance to FOM was mapped in a new BQ 
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population that included rfol, which is a loss-of-function al- 
lele of Atlg79670 resulting from a T-DNA insertion in and 
deletion of coding sequence in the Col-0 genetic back- 
ground [1,24]. The same crossing scheme that generated 
the original Col-0 and Ty-0 (C-T) BQ population, was 
used to generate the new rfol and Ty-0 (r-T) BQ popula- 
tion with the exception that rfol replaced wild type as the 
Col-0 parent: Crossing rfol and Ty-0 produced the Fi hy- 
brid that was then backcrossed to Ty-0. Differences in 
quantitative resistance in the new r-T and original C-T pop- 
ulations would include the contribution of Atlg79670. 

As in the C-T population, resistance to FOM segre- 
gated in the r-T population as a polygenic trait, and 
most BCi hybrids exhibited resistance that was inter- 
mediate to that of either parent [1]. Wilt disease in the 
Fx hybrid, Ty-0 parent and 190 BCi hybrids was evalu- 
ated using a health index (HI), an ordinal scale from 0 
(dead) to 5 (unaffected), described in Methods. At 18 days 
post infection (dpi), a broad distribution of HI scores reg- 
istered the breadth of disease resistance among BQ hy- 
brids and presumably the diversity of resistance genotypes 
(Figure lc). In contrast, the parents were consistently ei- 
ther resistant or susceptible. Most Fi hybrids (Figure la) 
as well as a minority of BCi hybrids (Figure lc) exhibited 
only mild symptoms (with a HI score > 3); and, at the 



opposite extreme, most of the Ty-0 parents (60 percent, 
Figure lb) as well as 10 percent of BCi hybrids were dead 
(Figure lc). Thus, segregation of resistance among BCi 
hybrids was inconsistent with monogenic inheritance as a 
single locus would have given a 1:1 segregation ratio in 
the backcross, i.e. one plant as resistant as the ¥ 1 hybrid to 
one plant as susceptible as the Ty-0 parent. 

Genome-wide linkage of 40 CHR markers 

To expedite the mapping of resistance, methodology to 
genotype BQ hybrids was developed with efficiency and 
economy in mind. Previously, RFO QTLs were mapped 
in the C-T population using the genome-wide genotype 
of 24 SSLP markers distributed over the five chromo- 
somes of A. thaliana [1]. However, genotyping one SSLP 
in one BCi hybrid from one PCR sample is a prohibitive 
bottleneck in analysis. For instance, if the same 24 SSLPs 
were used to genotype the 190 FOM-infected BCi hy- 
brids in the r-T population, the effort would entail pro- 
cessing no fewer than 4,560 PCR samples. Instead, as 
described in Methods, the genome-wide genotype of 40 
marker loci in each BCi hybrid was obtained from just 
three multiplex PCR samples. In comparison to genotyp- 
ing with SSLPs, the new approach gave genome-wide ge- 
notypes of BCi hybrids that were comprised of two-thirds 
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more markers and obtained in one-eighth as many PCR 
samples. 

The phenotype of the 40 CHR markers was dominant, 
and primer pairs for CHR markers directed PCR amplifi- 
cation of marker sequence from Col-0 DNA and not from 
Ty-0 DNA (Figure 2a). DNA products corresponding to 
as many as 14 markers were amplified in a single multi- 
plex PCR sample and then separated by size using stand- 
ard agarose gel electrophoresis, as shown for the three 
multiplex PCR samples of five representative BCi hybrids 
in Figure 2a. Because BCi hybrids were either Col-O/Ty-0 
(C/T) or Ty-O/Ty-0 (T/T) at any locus, genotypes were 
scored according to whether PCR-amplified marker DNA 
was present (C) or absent (T), respectively (Figure 2b). 



Genetic linkage between CHR markers in both r-T 
and C-T populations was consistent with the proximity 
and order of marker sequences in the Arabidopsis refer- 
ence genome (version TAIR10, www.arabidopsis.org). 
Genome-wide genetic maps corresponding to recom- 
bination frequencies in r-T and C-T populations are 
shown in Additional file 1: Figure SI and Additional file 2: 
Figure S2, respectively. In the r-T population, marker in- 
tervals had mean, median and total genome distances of 
15.8 centiMorgan (cM), 12.9 cM and 551 cM, respectively, 
while individual marker intervals ranged from 4.8 cM to 
27.1 cM. (See Additional file 3: Table SI for recombination 
frequencies and genetic distances of all intervals). In the 
original C-T population, 39 dominant markers and one 
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Figure 2 Genome-wide genotyping with CHR markers, (a) Multiplex PCR products for 40 Col-0-specific dominant markers were size-separated 
by agarose-gel electrophoresis and stained with ethidium bromide. Sizes in basepairs (bp) for the DNA ladder (leftmost lane) are at left. Marker 
DNA was PCR-amplified from the Col-0 and Ty-0 F, hybrid (C/T), accession Ty-0 (T/T) and C-T BC, hybrids, 5A2, 5A3, 5A4, 5A6 and 5A7. Lines to 
the right indicate the expected positions of marker bands. Markers are named CHRx.n, where x is the chromosome and n is the relative position 
on the chromosome, (b) Genotypes of markers in five BC, hybrids are from banding phenotypes in (a). Markers are ordered with respect to their 
position on chromosomes. Genotype C/T (C) is shown with black on white type, and genotype T^ (T) is white on black type. 
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SSLP marker (in place of the linked CHR2.4 marker) had 
mean, median and total genome distances of 14.1 cM, 
14.5 cM, and 516 cM, respectively (See Additional file 4: 
Table S2 for recombination frequencies and genetic dis- 
tances of all intervals). 

Reliability of CHR markers 

There was concern that dominant CHR markers would 
not be as reliable as codominant SSLP markers. The ab- 
sence of marker DNA, which is the phenotype of geno- 
type T/T, could be the false negative result of insufficient 
PCR amplification of Col-0 DNA from genotype C/T in 
which case a genotype of C/T would be miscalled as T/T. 
The codominant SSLPs, on the other hand, were safe- 
guarded from false negative miscalls because marker 
primers direct the amplification of Ty-0 DNA in all sam- 
ples, confirming that PCR was productive. 

Results with SSLPs and CHR markers were compared 
in the C-T population. As expected, half of genotypes at 
SSLP markers (50.4 percent with a standard deviation of 
3.4 percent) and half of genotypes at CHR markers (50.9 
percent with a standard deviation of 3.0 percent) were 
T/T, so neither codominant SSLPs nor the Col-O-specific 
CHR markers were prone to give an excess of T/T. 

The reliability of dominant markers was further scruti- 
nized by examining recombination in a dataset that 
combined the genotypes of 39 dominant markers and 24 
SSLP markers in the C-T population. Miscalled marker 
genotypes would exaggerate the number of instances of 
crossovers in adjacent marker intervals because tightly 
linked markers usually share the same genotype. The 
mean recombination frequency in intervals separating 
the 63 markers was 8.7 percent, so pairs of adjacent in- 
tervals were expected to have crossover events only once 
or twice among 234 BCi hybrids. A miscalled genotype 
would appear to be flanked, in most cases, by markers 
with opposite genotype and thus by intervals with spuri- 
ous crossovers. However, instead of having an excess of 
adjacent double crossovers, the combined marker dataset 
had a clear deficit of linked double crossovers (Figure 3a). 
A total of 80 double crossovers were predicted from the 
sum of the products of recombination frequencies in adja- 
cent intervals, whereas crossovers in adjacent intervals 
were observed in just 18 instances. Importantly, double 
crossovers flanked a similar proportion of dominant 
markers (10) and SSLPs (8). Thus, dominant markers were 
no more likely than SSLPs to have genotypes that were 
different from the genotypes of both flanking markers. In 
addition, the number of crossovers in two intervals was 
expected to decline as the number of marker intervals 
separating crossovers increased, whereas the observed 
number of double crossovers increased with separation of 
crossover events (Figure 3a). 
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Figure 3 Crossover interference, (a) Among 63 markers in 234 
FOM-infected C-T BCi hybrids, the number of expected (open bar) 
or observed (filled bar) crossovers is shown at the indicated distance 
between crossovers, (b) Coincidence of crossovers is the observed 
crossover frequency in two marker intervals divided by expected 
crossover frequency. Observed and expected frequencies are equiva- 
lent at 1 (dashed line). Crossovers are separated by number of 
markers (estimated distance in cM): 1 (9), 2 (18), 3 (26), 4 (35), 5 (44), 
6 (53) and 7 (62). The estimated distance between crossovers is the 
mean distance between adjacent markers (87 cM) times the number 
of markers. 
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Crossover interference 

The deficiency of linked crossovers was explained by 
crossover interference, which is observed in A. thaliana 
[25]. Coincidence, a measure of crossover interference, 
is defined as the observed frequency of crossovers in 
two marker intervals divided by the product of recom- 
bination frequencies in the same intervals [26]. Positive 
interference has a coincidence value of less than one and 
indicates that a crossover in one interval inhibits cross- 
over in the other interval. In Figure 3b, positive interfer- 
ence was observed when crossovers were separated by 
less than 36 cM. Coincidence values near one or greater 
than one indicate no interference or negative interfer- 
ence, respectively. A transitional negative interference, 
which is a common observation when positive interfer- 
ence is present, was apparent when crossovers were sep- 
arated by roughly 53 cM in Figure 3b [26]. Moreover, 
clear deficiencies of linked double crossovers were ob- 
served in all BCi populations examined here. Only 19 to 
33 percent of expected double crossovers in adjacent 
marker intervals were in fact observed in the five Arabi- 
dopsis chromosomes. (For chromosomal distribution of 
expected and observed double crossovers, see Additional 
file 5: Table S3.) 

No RFOl QTL without At1g79670 

In the r-T population, association of resistance at CHR 
markers was evaluated using the Mann-Whitney rank 
sum test as previously described in [1]. Briefly, BCi hy- 
brids in the r-T population were ranked, from most 
susceptible to most resistant, according to HI scores. At 
each marker, a standardized statistic Z enumerated the 
separation of ranks of BCi hybrids that were C/T and 
T/T, and the sign and magnitude of Z indicated the 
direction and strength of genetic association. Specifically, 



resistance was found to have significant correlation with 
genotype C/T when Z was greater than 3.3 and with T/T 
when Z was less than -3.3 (when p < 0.05, according to 
permutation tests). 

In the r-T population, no QTL with major effect was 
detected on chromosome 1, though both RFOl and 
RF02 are located on chromosome 1 and make substan- 
tial contribution to resistance in the C-T population 
(Figure 4) [1]. In the r-T population, the correlation of 
resistance with the Col-0 alleles of RFOl- and RF02- 
linked markers, respectively CHR1.9 and CHR1.3, lacked 
statistical significance (Figure 4). Thus, rfol abolished 
the major contributions of RFOl and RF02. In the C-T 
population, RF02's strong association with resistance 
among plants that are C/T at RFOl is absent among 
plants that were T/T (Figure 5a) [1]. In the r-T popula- 
tion, 7?F02-linked markers had insignificant association 
with resistance whether BQ hybrids were C/T or T/T at 
RFOl (Figure 5b). 

rfol also suppressed the apparent interactions between 
RFOl and either the (CHR5.6-linked) RF06 or (CHR4.2- 
linked) RF04. From prior work, resistance is associated 
with two loci on chromosome 5: RF05 gives resistance 
that is independent of RFOl while RF06 is only evident 
among BCi hybrids that are also C/T at RFOl (Figure 5a) 
[1]. In the r-T population, CHR5.6 lacked significant 
association with resistance among BCi hybrids with or 
without RFOl-C (Figure 5b). Similarly, an apparent inter- 
action between RFOl and (CHR4.2-linked) RF04 was not 
evident in the r-T population, whereas significant associ- 
ation of resistance at RF04 in the C-T population is evi- 
dent only among plants that also have RFOl-C (Figure 5a) 
[1]. In the r-T population, marker CHR4.2 was associated 
with a major QTL without regard to the genotype of 
RFOl (Figure 5b). 
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As previously observed, RF03 and RF05 expressed re- 
sistance that was independent of RFOl [1]. In fact, RF03 
and RF05 had stronger correlation with resistance in the 
r-T population than in the C-T population -compare 
peak Z values at 7?F03-linked (CHR3.3) and RF05- 
linked (CHR5.3) in Figure 4. Excluding RF03, RF04 and 
RF05, there was no other significant association with re- 
sistance to FOM. 

RFO QTLs are pathogen-specific 

To examine whether the same or different RFO loci pro- 
vided resistance to different F. oxysporum pathogens, re- 
sistance was investigated in a third BCi population that 
was instead infected with FOCI. The HI scores of 200 
FOCl-infected BQ hybrids and the two parental acces- 
sions, Col-0 and Ty-0, at 16 dpi are shown in Figure Id, e 
and f, respectively. Both parental accessions exhibited par- 
tial resistance to FOCI and had median HI scores of 3.5 
(Figure Id, e). BCi hybrids exhibited a broader range of 
symptom severity than their parents (Figure If): 17 per- 
cent of BCi hybrids were unaffected (HI = 5.0) while all 
parents exhibited at least mild symptoms; and, 15 percent 
of BCi hybrids exhibit more severe symptoms that either 
parent (HI < 2.0). Thus, a third of FOCl-infected BQ hy- 
brids expressed an extreme phenotype that was not seen 
in either parent. 



RF07 confers resistance to FOCI 

A genome-wide genetic map derived from the recombin- 
ation frequencies between CHR markers in the FOCl- 
infected C-T population was consistent with the order of 
marker sequences in the TAIR10 reference genome (See 
Additional file 6: Figure S3 for the genome-wide genetic 
map). Intervals between markers ranged from 4.0 to 24.9 
cM, and mean, median and total genome distances were 
13.5, 12.8 and 472 cM, respectively (See Additional file 
7: Table S4 for recombination frequencies and genetic 
distances of all marker intervals). 

Association with resistance to FOCI was evaluated at 
40 CHR markers. For the sake of comparison, Z statis- 
tics at markers in FOCl-infected and FOM-infected C-T 
populations are juxtaposed in Figure 6. A single major 
effect QTL at marker CHR5.7 (Z = -8.77) associated geno- 
type T/T with strong resistance to FOCI. Because all pre- 
vious RFO QTLs correlated resistance with genotype C/T 
and CHR5.7 was not previously associated with resistance, 
this QTL was new and was named RF07. Among F 2 off- 
spring of Col-0 and Ty-0, genotype C/C at the RF07- 
linked SSLP CIW9 was more susceptible to FOCI than 
genotype C/T, indicating that Col-0 and Ty-0 alleles of 
RF07 express incomplete dominance (Figure 7). 

Previously, RFOl was shown to confer resistance to 
FOCI as well as FOM [1]. In the FOCl-infected C-T 
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Figure 6 Association of resistance to FOCI. In FOM-infected (circles) and FOCI-infected (triangles) C-T populations, test statistic Z enumerates 
the correlation of resistance and marker genotype, and lines connect linked values. See Figure 4 for description of plot details. 



population, resistance associated with 7?FOi -linked CHR1.9 
had questionable significance (Z = 2.54, p = 0.28). How- 
ever, among BCi hybrids that were heterozygotes (C/T) 
at CHR5.7, which minimized the contribution of RF07, 
the association of resistance with RFOl was significant 
(Figure 8). 

Discussion 

It was conceivable that more than one gene might be 
responsible for different aspects of the RFOl QTL. How- 
ever, QTL analysis that included rfol was consistent with 
the simplest explanation: A single gene was responsible 
for the major effect of RFOl and also for apparent inter- 
actions with RF02, RF04 and RF06. RF02 and RF06 
were undetected in the r-T population while the resist- 
ance of RF04, which was dependent on RFOl-C in the 
C-T population, was independent of RFOl in the r-T 
population. 

Why RF04, which only attained significance with 
RFOl-C in the C-T population, was a major QTL in the 
r-T population lacking RFOl-C is difficult to explain. 
Possibly, the expression of RF04 was influenced by sub- 
tle differences in the progression of wilt disease or envir- 
onmental factors as the C-T and r-T populations were 
similarly infected on separate occasions. Also, the par- 
ents, which were nominally from the same Col-0 and 
Ty-0 accessions, might have been genetically (or epige- 
netically) dissimilar as separate crosses generated the 
two populations. 

Overall, results obtained from independent FOM- 
infected populations were consistent, and QTLs in the 
new r-T population were coincident with the previously 
detected RF03, RF04 and RF05 in the C-T population. 
In fact, the association of resistance at the three QTLs 
appeared stronger in the r-T population. Col-0, which was 
the source of resistance, rarely exhibits wilt symptoms 
when infected with FOM and presumably expresses more 
than sufficient resistance. As symptom severity is difficult 



to discriminate among the more resistant plants, loss of 
RFOl-C in the r-T population undoubtedly improved the 
evaluation of disease in BCi hybrids and thus the detec- 
tion of RFO QTLs. 

The map position and source of QTLs detected in 
FOCI- and FOM-infected populations suggest that 
quantitative resistance to F. oxysporum is predominantly 
specific to the infecting forma specialis. Remarkably, re- 
sistance to FOCI was strongly associated with a single 
new QTL, RF07, though Col-0 and Ty-0 expressed simi- 
lar partial resistance to FOCI. In addition to RF07, a re- 
gion on chromosome 1 had a marginal association with 
resistance, and RFOl -linked markers did attain signifi- 
cant association when BQ hybrids were heterozygous 
for RF07. However, previous work clearly shows that 
rfol and transgenic RFOl affect resistance to both FOCI 
and FOM. Thus, while RFOl may have a non-specific 
role in resistance to the three crucifer-infecting formae 
speciales, it appears that RFOl also has a much stronger, 
specific effect on resistance to FOM. 

Only dominant traits from the donor parent are 
expressed in a BCi population, so RFO alleles of Col-0 
that were recessive to alleles of Ty-0 would not be de- 
tected. Nevertheless, the strong resistance of Fi hybrids, 
Col-O/Ty-0 and rfol /Ty-0, suggests that resistance to F. 
oxysporum is in large part a dominant trait. Resistance 
associated with RF07 was confirmed in F 2 progeny, and 
positional cloning has identified single genes that are re- 
sponsible for three RFO QTLs, RFOl, RFOl and RF03 
[1,23,27]. Of the four confirmed RFO loci, three QTLs, 
RFOl, RF03 and RF07, express incomplete dominance. 

Nucleotide sequences of both resistant and susceptible 
alleles of RFOl, RF02 and RF03 encode apparently func- 
tional, full-length membrane-spanning receptor proteins. 
Thus, competition between or interference by the products 
of the two alleles, rather than gene dosage of the resistant 
Col-0 allele, might explain the incomplete dominance of 
natural RFO alleles. Because physical interactions between 
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(HI score of 0 or 1, black-filled bar) or had intermediate resistance (H 
scores of 2 or 3, gray-filled bar). M-W test indicates that symptom 
severity in C/C and C/T (p = 0.005) or in C/C and T/T (p = 0.0006) 
was dissimilar. 



RLKs and a RLP are critical for signaling in plants [28], 
genetic interaction between RFOl and RF02 might be evi- 
dence for the direct interaction of the corresponding RFOl 
RLK and RF02 RLP. However, because resistance is a 
complex phenotype, involving processes that occur at dif- 
ferent sites in the host and at different times in the infec- 
tion cycle, the observed genetic interaction might reveal 
the priority of RFOl before RF02 without direct inter- 
action. In addition, when neither resistant nor susceptible 
allele is the null allele, interpretation of genetic interaction 
is ambiguous. For example, results with rfol clearly impli- 
cate RFOl-C in the RF01-RF02 interaction, however, it 



remains unclear whether RFOl-C suppresses resistance of 
RF02-T or enhances resistance of RF02-C. 

Routine QTL analysis in A. thaliana is limited to natural 
traits that distinguish parents of existing RI populations; 
otherwise, the generation of new RI lines represents a sub- 
stantial investment of time and effort [4,5]. In the mean- 
time, there is an increasing availability of whole genome 
sequence, from the 1001 Genomes Project for example, 
that makes the sequence diversity in hundreds of Arabi- 
dopsis accessions accessible. As shown here, the mapping 
of traits that distinguish any two sequenced accessions, in- 
cluding mutant genotypes, can be conceive and complete 
in six months using BCi populations. 

Mapping in BCi populations can be a routine procedure 
when genotyping is efficient, accessible and economical. 
In this regard, available whole genome sequence from the 
1001 Genomes Project was an invaluable resource for 
identifying accession-specific polymorphisms [29]. Primer 
sequences for Col-0-specific dominant markers were read- 
ily selected from genome sequence reported to be poly- 
morphic in Col-0 and Ty-0. In the same way, dominant 
markers could be designed to distinguish any two se- 
quenced accessions. In fact, we have reused most of the 
Col-0-specific markers for genotyping BCi populations 
from crosses between Col-0 and accessions Zdr-1 or 
Kondara (unpublished observation). 

The methodology for genotyping was designed with effi- 
ciency and economy in mind. Starting from crude leaf 
preparations, multiplex PCR DNA of 40 dominant markers 
was amplified by just three sets of multiplex PCR primers 
and visualized using standard agarose-gel electrophoresis. 
The 200 FOCl-infected BQ hybrids were genotyped 
genome-wide with little more than six 96-well plates of 
PCR samples. Markers in a multiplex PCR sample appeared 
as a ladder of bands in agarose gels when all markers were 
present. Because annealing of marker primers distinguished 
the Col-0 and Ty-0 genotypes, markers could be arbitrarily 
assigned sequence lengths that appeared as regularly spaced 
bands in agarose gels. 

Results obtained with dominant markers were as reli- 
able as results from codominant SSLP markers [1]. No 
unforeseen PCR products were amplified when as many 
as 14 primer pairs were combined in multiplex PCR, 
and no primer pairs that were confirmed singly subse- 
quently failed when combined with other primer pairs. 

In theory, an RI population has roughly twice as many 
crossovers as a BCi population [30]. However, the add- 
itional recombination in RI lines remains largely un- 
appreciated unless a high density of DNA markers are 
used to genotype RI lines [31]. During the inbreeding cy- 
cles that generate RI lines, crossovers tend to accumu- 
late at linked sites, and thus recombination in RI lines 
has the appearance of negative interference. High- 
resolution analysis of breakpoints in 98 Col-0/Ler RI 
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lines found, for example, that 17 percent of intervals be- 
tween crossovers contained just one gene [32]. 

For genome-wide linkage analysis in BCi populations, 
40 evenly spaced markers should be sufficient to capture 
most recombination. As already mentioned, just one set 
of homologous chromosomes in BCi hybrids is recom- 
binant. With an average marker separation of 15 cM, I 
estimate that just seven percent of crossovers went un- 
detected in the BCi populations because approximately 
three quarters of expected double crossovers in marker 
intervals would be suppressed by positive interference. 
In addition, I took advantage of significantly higher re- 
combination in male meiosis (as compared to female 
meiosis) in A. thaliana when generating the BCi hybrids 
[33]: Ty-0 was the female parent in the backcross while 
the Fi hybrid, which was the source of recombinant 
chromosomes, was the male parent. 

Number of crossovers, or amount of recombination, 
has little bearing on whether a lone QTL is detected 
[4,13]. Rather, recombination frequency affects the reso- 
lution of map position of a QTL, and less recombination 
would more poorly resolve multiple QTLs in proximity 
on a chromosome. The detection of two or more linked 
loci could be suppressed if the loci that remain unre- 
solved express opposing effects on a trait. Indeed, an 
example of two opposing QTLs for growth rate within 
an interval of 210 kbp has been reported in Col-0/Ler 
recombinants [34]. 

QTL mapping is just the first step in the identification 
and characterization of the genes underlying traits. In 
this regard, mapping in BCi populations is also advanta- 
geous because individual (or specific combinations of) 
QTLs can immediately be reevaluated and fine-mapped 
in progeny of selfed BQ hybrids. Even after a potentially 
lethal test, such as resistance to FOM, I was able to collect 
seeds from 144 of 236 tested C-T BCi hybrids. Although 
half of the genome in BCi hybrids was heterozygous, on 
average, seeds were collected from 16 BQ hybrids that 



were largely homozygous Ty-0 and heterozygous in just 
four or fewer chromosomal intervals representing 30 per- 
cent or less of the genome. RFO QTLs in these heterozy- 
gous intervals would again segregate in progeny. 

Conclusions 

Genome-wide mapping of quantitative Fusarium wilt re- 
sistance was expeditious and reproducible in BCi recom- 
binant populations of A. thaliana. In two independent 
BCi populations, resistance to FOM was associated with 
QTLs RF03, RF04 and RF05. Because the resistance of 
RFOl, RF02 and RF06 was absent in the BQ population 
that included rfol, the major effect and epistatic interac- 
tions of RFOl were solely attributed to Atlg79670, the 
gene disrupted in rfol. In a third BQ population, resist- 
ance to a second pathogen FOCI was instead associated 
with RF07, a new major effect QTL. Pathogen-specific 
RFO QTLs were largely responsible for resistance to the 
two pathogens, FOM and FOCI. 

Methods 

Growing A. thaliana 

Seeds of Ty-0 (CS6768) and rfol (Salk_077975) were ob- 
tained from the Arabidopsis Biological Resource Center. 
Seeds were surface-sterilized in 10% household bleach 
and 0.1% Triton X-100 for 15 min, rinsed 3 times in 
sterile water. Seeds were sown on peat pellets (Jiffy-730, 
Grower's Solution Inc., Cookeville, TN) or first germi- 
nated on plant nutrient agar (PNA) before transplanting 
[1]. Plants were arrayed in flats (1' x 2') 5 rows by 10 
columns and designated: first by flat (1 through 6), sec- 
ond by row (A to E) and third by column (1 through 
10). Plants were grown under medium intensity cool 
white fluorescent lighting (100 to 140 umoles rrf 2 sec' 1 ) 
for a 12-hr daylength at 25 to 28°C and irrigated with 
water or fertilizer (PlugCarePlus, Greencare Fertilizers, 
Inc., Kankakee, IL). 
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Infection with F. oxysporum 

Fusarium oxysporum forma specialis conglutinans race 1 
(FOCI, isolate 777) and Fusarium oxysporum forma 
specialis matthioli (FOM, isolate 726) are from P.H. 
Williams by way of H.C. Kistler [3,35]. F. oxysporum cul- 
tures were stored at -80°C in 50% glycerol, grown on 
Czapek Dox medium (Oxoid Ltd., Hampshire, England) 
and harvested as described in [1]. Starting with an excess 
of 3-week old plants, 200 C-T BQ, 25 Col-0 and 25 Ty- 
0 plants with comparable sizes were infected with FOCI; 
and, 190 r-T BCi, 25 rfol, and 15 Ty-0 plants were in- 
fected with FOM. Plants were irrigated with an excess of 
washed conidia (2 x 10 6 conidia ml/ 1 ). The FOM-infected 
r-T population was scored 1 1 days post infection (dpi) for 
three early symptoms: (i) stunting of leaves, (ii) leaf epi- 
nasty and (iii) anthocyanin accumulation, using a gradu- 
ated scale of 1 (severe) to 4 (unaffected). At 18 and 23 dpi, 
infected plants were scored using a health index (HI), 
which is the same as the disease index (DI) in [1], ranging 
from 0 (dead plants) to 5 (unaffected plants) in intervals 
of 0.5. The FOCl-infected C-T population was similarly 
scored on 10, 13 and 16 dpi. At the final time point, plants 
were rank ordered: For the FOCl-infected C-T popula- 
tion, each flat of 40 plants was ranked separately, from 1 
(most susceptible) to 40 (most resistant); and, for the 
FOM-infected r-T BC population, all 190 plants were 
ranked together, from 1 (most susceptible) to 190 (most 
resistant). Infection and scoring of the FOM-infected C-T 
BC population is described in [1]. (See Additional file 8: 
for spreadsheets with phenotypic data of all three BCi 
populations). 

Genotyping with CHR markers 

In proportion to physical and genetic lengths of chromo- 
somes in The Arabidopsis Information Resource (TAIR, 
www.arabidopsis.org), 10, 6, 8, 7 and 9 CHR markers were 
distributed on chromosomes 1, 2, 3, 4 and 5, respectively. 
On each chromosome, two markers were placed close to 
the telomeres, and nucleotide positions for remaining 
markers were spaced at regular intervals in the reference 
Col-0 sequence of TAIR10. 

At the approximate nucleotide positions of markers, 
marker sequences were reference sequences that were 
classified as highly diverged, or "unsequenced", in whole 
genome sequencing of Ty-0 (http://signal.salk.edu/ 
atglOOl). Appropriate pairs of primer sequences were 
selected in the highly diverged reference sequence using 
Primer3Plus software, according to recommendations 
of the QIAGEN Multiplex PCR Handbook (Qiagen Inc., 
Valencia, CA) [36]. DNA products of 13, 13 or 14 
markers were simultaneously amplified by multiplex 
PCR using three sets of PCR primers. Each set of multi- 
plex PCR primers were designed to give a logarithmic 
progression of DNA product sizes, ranging from 200 bp 



to 650 bp in length, which gave regular spacing of 
marker bands when products were size-separated by 
agarose gel electrophoresis. Primer sequences and gen- 
omic locations of PCR primers are in Additional file 9: 
Table S5. Sizes and order of DNA products for each set 
of primers are in Additional file 10: Table S6. PCR amp- 
lification was performed using the QIAGEN Multiplex 
PCR kit according to the protocol for microsatellite 
loci. A reaction volume of 5 uL included 1 uL crude leaf 
DNA preparation, 2.5 uL 2x QIAGEN Master Mix, 1 uL 
lOx primer mix (containing 2 uM of each oligonucleotide 
primer), and 1.5 uL water. Amplified PCR products were 
separated by gel electrophoresis in 2% agarose. Crude leaf 
DNA preparations were prepared according to [37]. See 
Additional file 8: for spreadsheets with genotypic data for 
all markers in all BCi populations. Genotypic data for sim- 
ple sequence length polymorphism (SSLP) C4H, from a 
prior study [1], replaced CHR2.4s in the analysis of the 
FOM-infected C-T population. 

Genetic distances between markers were calculated 
using the Kosambi mapping function [13]. Genetic 
linkage supported the presumed physical linkage of 
markers in the three BCi populations. Linkage data for 
markers in the three mapping populations are provided 
in Additional file 3: Table SI, Additional file 4: Table 
S2 and Additional file 7: Table S4. 

Testing association of wilt resistance 

A BCi population of n plants was rank ordered accord- 
ing to HI scores, from 1 (most susceptible) to n (most 
resistant). Ranking gave priority to later HI scores over 
earlier HI scores. Rank distributions of the two possible 
genotypes C/T and T/T were compared using the Mann- 
Whitney (M-W) test, the results of which were expressed 
as a standardized statistic (Z), the standard deviation units 
separating the mean ranks of the two genotypes. (See 
Additional file 11: Table S7 with values of Z at CHR 
markers in the three BCi populations.) For a major effect 
QTL, threshold values of Z for the three BCi populations 
were determined by permutation tests with 10,000 trials 
[18]. From the distribution of highest Z values in trials, 
the threshold value of Z at p = 0.01 (Z 0 0 i) was 3.86; and, 
Zn.o5 was 3.36; and, Z 02 o was 2.80. In r-T population, 
Zn.oi = 3.73; and, Z 0 .os = 3.27; and, Z 0 .2o = 2.77. In FOCl- 
infected C-T population, Z 00l = 3.63; and, Z 005 = 3.16; 
and, Z 0 .2o = 2.68. Probability threshold values of Z were 
also determined for QTLs conditioned by genotype at a 
major QTL [20]. For FOCl-infected C-T population, the 
M-W test was performed on subpopulations that were ei- 
ther genotype T/T (Z 0 .ni = 3.54; and, Z 0 .os = 3.10; and, 
Zn.20 = 2.67) or C/T {Z 001 = 3.51; and, Z 005 = 3.08; and, 
Z0.20 = 2.63) at 7?F07-linked CHR5.7. For FOM-infected 
plants, subpopulations were tested that were either T/T 
(for r-T population, Z 001 = 3.56; and, Z 005 = 3.13; and, 
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Zo.20 = 2.67; and, for C-T population, Z 0 .oi = 3.55; and, 
Zn.os = 3.12; and, Z0.20 = 2.66) or C/T (for r-T population, 
Zom = 3-52; and, Z 0 05 = 3.12; and, Z 0 . 2 o = 2.64; and, for C- 
T population, Z 0 .oi = 3.88; and, Zn.os = 3.34; and, Z 0 .2o = 
2.81) at RFOl-lmked CHR1.9. 
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